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Abstract 



We derive a trace formula for Y^n^nnSnn- ■ ■^{E — where A 
is the diagonal matrix element of the operator A in the energy basis 
of a chaotic system. The result takes the form of a smooth term plus 
periodic-orbit corrections; each orbit is weighted by the usual Gutzwiller 
factor times ApBp . . ., where Ap is the average of the classical observable 
A along the periodic orbit p. This structure for the orbit corrections 
was previously proposed by Main and Wunner on the basis of numerical 
evidence. 
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I. INTRODUCTION 



In a recent paper [|l|, Main and Wunner introduced the weighted density of states 

p"^^'""'-^ =Y.A^An - ■ - SiE - E„) . (1) 

n 

Here A, i?, . . . are operators with smooth classical limits (Weyl symbols), and Ann = {n\A\n) 
is the diagonal matrix element of A in the energy basis. This is a simple generalization of 
p(A) ^ J2nAnnS{E — En), which has been studied extensively pj-Q. For chaotic systems, 
Main and Wunner proposed that 

P^^'^'-^ = pr'-^ + ^E^.5p---^., (2) 

TTll p 

where the sum is over all primitive periodic orbits p with energy E, and 

Aj,= - rA{-Xp{t))dt (3) 



Tp Jo 

is the average value of the Weyl symbol of A along the orbit; here Tp is the period of the 
orbit, and X = (q, p) denotes both coordinates and momenta. Also, 

oo ^ ^i{Sp/h-fj.pTT/2)r 

"'"'''Sldet(Af;-/)|./. W 

is the Gutzwiller weight factor; Sp, fip, and Mp are respectively the action, Maslov index, and 
monodromy matrix of the orbit. The first term on the right-hand side of Eq. represents 
the part that remains smooth in the semi classical limit; it should be 0{h~^), where / is the 
number of freedoms. However, Main and Wunner do not give an explicit formula for it. 

If we set -B = . . . = J in Eq. ([l|), and hence Bp = ... = 1 in Eq. (H), then we recover the 
trace formula for p^'^^ If we set A = I as well, we recover the original Gutzwiller trace 

formula for the density of states . This is the essential motivation of Main and Wunner 
for Eq. (H). They provide strong numerical evidence in favor of it, but they do not give an 
analytic derivation. 

In this paper, we remedy this situation by deriving Eq. (|]) from a generalization of a 
trace formula originally due to Wilkinson |^ (see also Furthermore we provide an 

explicit expression for the smooth term. 



II. ANALYSIS 

We first consider the case of two operators, A and B, and extend the results to an 
arbitrary number in Section Following Wilkinson , we define 

S{E, A) = J2 AnmBmn6l{E - \{En + E.^))62{^ " (K " ^m)) • (5) 

nm 

Here 5i{E) and 52{E) are smeared delta-functions. Rigorous results concerning S{E,A) 
have been proven in the case that the Fourier transforms of these smeared delta functions 
have compact support @J^. We will therefore make the simple choice 
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where Tj, i = 1,2, is a time cutoff. Our results will come from various manipulations of 
S{E, A) with A = 0. 
We begin by writing 

S{E, 0) = 52(0) E AnnBnJi{En - E) 

n 

{E-\{Er, + E^)) , (8) 

where Unm = (-E'n — Em)/h. The first term on the right-hand side is the one we want; except 
for the factor of 52(0) = it is the same as the right-hand side of Eq. (0), in the limit 

as Ti ^ 00. To get rid of the unwanted second term, we take T2 to be much greater than 
the Heisenberg time = 2TThpQ] here 

is the Weyl formula for the mean density of states. If T2 ^ then we typically have 
||^nm|T2 ^ 1. In this case, sin(cij„mr2) varies erratically as n and m are varied. Furthermore 
the factor of 1/ {En — E^) can be written as pQ/{n — m) = [271% /th) /(n — m), up to a factor 
which also varies erratically. We then have 

S{E, 0) = 4 ^ AnnBnJi{En - E) 

+ 2?^ ^ ^";'^_-j^- 5,(i?-i(i?„ + i?^)), (10) 

where we can think of Rnm as a random number. Provided that \Amn\ and \Bmn\ do not 
tend to increase as \m — n\ increases (in general a decrease is to be expected), the sum in 
the second term should quickly converge. Then we have 

— SiE, 0) = Y. AnnBnJliEn - E) + 0{th/t2) . (11) 
''2 n 

The first term on the right-hand side is the same as the right-hand side of Eq. (^, provided 
'Ti^ thi and the second term is small if T2 ^ th- 

We now wish to evaluate S{E, A) semiclassically. We first use Eq. (|^) in Eq. (|^) to get 

where we have defined 
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The key point is that we can write t') as a single trace, 

F{t, t') = Tr U{-t + \t')AU{t + \t')B , 



(14) 



where U{t) = e~*^*/^ is the time-evolution operator. 

To simplify our exposition, we temporarily make the (otherwise unnecessary) assumption 
that the Weyl symbols of A and B are functions of only the coordinates q and not the 
momenta p. We can then evaluate the trace by inserting two complete sets of position 
eigenstates, leading to 

F{t,t') = J dfq,dfq2{qi\Ui-t + lt')\q2) A{q,) {q^lUit + lt')\q^) B{q,) . (15) 
We now make use of the semiclassical approximation 0,|1O[ to get 

'd^q2 (q3|f/(-t + |t')|q2) ^(q2) (q2|f/(t + |t')|qi) 



E ^Path(q3, qi; t') A(qpath(t + It')) . 

paths 



(16) 



Here the sum is over all classical paths that go from qi at time zero to qs at time t', 
qpath(T) is the position reached at time r along a particular path, and i^^path(q3, qi; t') is the 
contribution of that path to the propagator (q3|f/(t')|qi) in the semiclassical limit. 

We now perform the integrals over d^qi in Eq. (|T^) and over dt' in Eq. (|T2[) by stationary 
phase p|-0. We get a contribution from zero- length paths (for which t' = at the point of 
stationary phase), and a sum over contributions from periodic orbits (for which t' = Tp at 
the point of stationary phase). The result is 



S{E,A) 



+^2 dt 



2nh 



-iAt/h 



PoCoit) + ^ E WpCpit) 



(17) 



where the sum is over all primitive periodic orbits with period less than ti. Also, we have 
introduced the energy-surface correlation function 



Co{t) 



1 



PoJ {2nh)f 
and the orbit correlation function 



6{E - H{X)) A{X{t))B{X) 



Cp{t) = - r dr A{Xp{r + t + lTp))B{Xp{r)) . 
r„ Jo 



(18) 



(19) 



Next, we must separate out a possible constant term in Co{t). To do so, we take the 
microcanonical average of ^(X) on a surface of constant energy E, 



Ao^- 



d'^fx 



PoJ {27rh)f 



S{E-H{X))A{X), 



(20) 



and define A{X) = A{X)-Ao and B{X) = B{X) - Bq. We then have Co{t) = AoBo + Co{t), 
where 
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Coit) 



6{E - H{X)) A{X{t))B{X) . 



Since the system is chaotic (and hence mixing), Co{t) ^ as t ^ ±00. We now have 



(21) 



SiE,A) = poAoBoS2{A) + 



+^2 dt 



-iAt/h 



1 



-T2 27!'h 

Next, we use the fact that Cp{t) is periodic in t with period Tp, which allows us to write 



(22) 



+00 

fc=— 00 



where 



We note in particular that 



pk 



Tp Jo 



dte-'^^'^''^^ Cp{t) 



pO 



ApBp . 



(23) 



(24) 



(25) 



Now using Eq. ( |23D in Eq. (|22D, we get 



^(E,A) = po^o5o52(A)+po 



+^2 dt 



2nh 



+ — ^ IWp ^pkS2{^-27lhk/Tp) . 

rp<Ti k=-oo 



(26) 



This result is a slight generalization of Wilkinson's P]. 

We now set A = 0. For the k = term in the sum over orbit modes, we have a factor 
of ^2(0) = T2/7!'h; for the A; 7^ terms, we use Eq. (|^). After dividing through by ^2(0), and 
using po = Tn/^irh, the result is 



^S{E,0) = poAoBo + ^ 

T2 2T2 



+T2 dt ~ 1 

KT^ C'o(t) + — E ^p^po 
T2 2Txn Txri 



Tp<Tl '2 fe^O 



Tp<n 

Tpk sin(27rA;r2/rp) 
27rA; 



(27) 
and 



If we assume T2 ^ Ti, the last term can be neglected. Comparing Eqs. (|TT|) and (^ 
recalling Eq. (^), we verify Eq. (|]). This is our key result. 

There is, however, an important caveat. We have taken both ti and T2 to be much 
greater than th ~ The rigorous treatment of on the other hand, requires xi 

and T2 to remain fixed as ^ — > 0, This suggests that a compromise of ti^2 ~ th might be 
optimal, which is consistent with other analyses of the trace formula [11-14]. To further 
investigate this issue, we consider the special case A = B. The magnitude of fluctuations in 
the values of the diagonal matrix elements of an operator A have been previously evaluated 
@J|J^,||,|1|], with the resuh that 
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pi^'^^=PoA^ + ^r"^C?o(t). (28) 

Here g = 2 ii the system is time reversal invariant, and = 1 if it is not. Eq. ( ^8]) is 
consistent with Eqs. ([l^) and (p7D if we set T2 = th/^q. We therefore conclude that, in the 
more general case where B, 

TT^ Coit) . (29) 

If we also set ti = r2, then the final term in Eq. (^) will be negligible for short periodic orbits 
(although it may become significant as the orbit period Tp approaches T2). This provides a 
theoretical explanation for the numerical results of Main and Wunner 

Eq. (|29| ) is an explicit formula for the smooth term. However, for ri^2 ~ th, it is an 
uncontrolled approximation, since the neglected terms in Eqs. ( p^lD and ( p7D are not formally 
suppressed. Nevertheless, the agreement of Eq. (pSD with the results of P,|6|j7| JT5|JT6[] leads 
us to believe that Eq. (P5|) is correct. 



III. EXTENSIONS 

We now return to Eq. (|I]) and consider a string of diagonal matrix elements of opera- 
tors, A, B, . . . , Z . We can find a trace formula for Eq. (0) by starting from 

F{h, . . . , t^) = Tr U{U)AU{h)B . . . U{tN)Z . (30) 

We then make Fourier transforms with respect to suitable linear combinations of the tj's to 
construct 

S{E, A2, . . . , A^r) = ^ AnmBml ■ ■ ■ h{.E — ^{En + Em + •••)) 



nml.. 



52(A2 - (K - Em))S3iAs - {Em - E{)) .... (31) 



A straightforward generalization of the analysis in Section then leads to a result analogous 
to Eqs. (0) and dg^), thus verifying Eq. (|). 

The leading contribution to the smooth term is of the form PqAq . . . Zq. There are also 
subleading contributions that depend on various energy- surface correlation functions. For 
example, for = 3 we have 

p(^'^'^) = poAoBoZo + ^ /^J dt [AoC^^it) + BoCl\t) + Z,C^^{t) 
I f g \^ f+TH f+TH ^^^^ 



+ I dUj dhC^^'{h,h). (32) 



Here CQ^{t) is given by Eq. 

C^^'ihM) = -j S{E - H{X)) A{X{t^))B{Xit,))Z{X) , (33) 

and again g = 2 if the system is time-reversal invariant, and (7 = 1 if it is not. This follows 
from requiring Eq. (|32|) to reproduce Eq. (^91) when we set Z = I. 
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IV. CONCLUSIONS 



We have formulated and verified a precise version of Eq. (Q), which was originally pro- 
posed by Main and Wunner [|l| for chaotic systems. Our derivation provides an analytic 
explanation for their numerical results. 

Main and Wunner also proposed an equation analogous to Eq. (|]) for integrable systems 



that generalizes the Berry- Tabor trace formula |T^. We have not attempted to derive this 



version, but clearly it would be of interest to do so. 
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